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摘 要 : 瞬时 射频 干扰 (RF1) 信号 是 最 难以 识别 并 消除 的 一 类 干扰 。 针 对 射电 天 文 观 
测 过 程 中 遇 到 瞬时 RF1 问 题 ， 本文 提 出 了 一 种 基于 RF1 信 号 峰值 的 特征 提取 方法 , 利用 该 方法 
生成 了 RF1 信 号 的 特征 模板 ， 并 基于 动态 时 间 规 整 (DTW) 算法 和 打分 策略 设计 了 特征 识别 算 
法 ， 实 现 了 RF1 信 号 识别 以 及 分 类 功能 。 利 用 ARRL 官 网 提供 的 RF1 数 据 完成 了 特征 模板 的 交叉 
验证 对 比 ， 测 试 结果 显示 RF1 特 征 模板 与 同类 别 RF1 相 似 度 高 ， 表 明 该 算法 能 够 对 RF1 进 行 识 
别 和 分 类 ， 该 算法 为 瞬时 RF1 的 特征 提取 以 及 识别 提供 了 一 种 新 思路 。 
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随 着 射电 天 文 技术 的 进步 ， 射 电 望 远 镜 向 大 口径 及 阵列 方向 发 展 。 中 国 500 米 口径 球面 
射电 望远镜 (Five-hundred-meter Aperture Spherical radio Telescope, FAST) 中 、 蓓 兰 低频 阵 
列 (the Low-Frequency Array, LOFAR) 外 以 及 美国 新 墨西哥 州 的 甚大 天 线 阵 (Very Large 
Array, VLA) 口 等 观测 设备 具备 高 灵敏 度 和 分 辨 能 力 。 在 观测 宇宙 中 微弱 天 文 信号 的 同时 ， 
干扰 信号 也 被 接收 ， 这 种 强度 比 微弱 天 体 信号 强 几 个 数量 级 的 人 为 射电 干扰 信号 (RFI) ， 
在 射电 天 文 数据 处 理 过 程 中 必须 得 到 有 效 抑制 。 

RFI 来 源 多 种 多 样 , 文献 [和] 将 RFI 分 为 外 部 RFI 和 内 部 RFI, 外 部 RFI 特 指 天 文 台 址 外 产生 
的 无 线 干扰 ， 常 见 的 有 人 造 卫星 外 (北斗 和 GPS 导 航 卫 星 等 )、 飞 机 、 附 近 的 基站 信和 号 以 及 广 
播 信号 等 ， 内 部 RFI 是 指 在 天 文 台 址 内 使 用 的 电子 设备 ， 如 各 类 数字 终端 设备 、 计 算 机 、 微 
波 炉 、 咖 啡 机 、 荧 光 灯 以 及 其 他 电子 仪器 等 产生 的 RFTI。 

射电 天 文学 发 展 过程 中 ，REFI 抑 制 一直 是 研究 热点 。 文 献 [6] 提 到 天 文 信号 通常 为 宽频 带 
且 往 往 在 时 间 尺 度 上 平滑 变化 ， 绝 大 部 分 RFI 与 天 文 信号 相 比 有 明显 的 差异 。 针 对 RFI 抑 制 ， 
文献 [7] 提 出 了 CUSUM 方 法 ， 该 方法 简单 、 高 效 ， 但 依赖 于 单 采 样 点， 对 噪声 比较 敏感 ; 为 
了 克服 这 一 局 限 性 ， 文 献 [8] 提 出 了 基于 计算 相 邻 采样 点 的 组 合 病 值 表示 算法 
—SumThreshold; 文献 [9] 利 用 奇异 值 分 解 来 解决 GMRT (ERRE: Giant Metrewave Radio 
Telescope 〉 宽 带 RFI 问 题 。 这 些 算法 本 质 上 基于 阔 值 实现 ， 处 理 幅 值 远 大 于 天 文 信号 的 RFI 
言 号 时 效果 较 好 ， 但 对 于 特征 与 天 文 信号 接近 的 RFI 信 号 ， 闵 值 法 可 能 会 误 删 天 文 信号 ， 导 
致 信息 丢失 ， 影 响 天 文 研究 的 准确 性 和 有 效 性 。 

瞬时 RFI 是 某 些 设备 或 系统 正常 运行 时 释放 的 射频 信号 。 继 电器 、 电 动机 和 开关 电源 都 
会 产生 瞬时 干扰 。 由 于 瞬时 干扰 信号 是 宽频 的 、 脉 冲 的 ， 而 且 持 续 时 间 非 常 短 ， 难 以 识别 来 
源 。 文 献 [10] 提 出 了 基于 字典 的 瞬时 RFI 分 类 算法 , 使 用 隐 马 尔 科 夫 模 型 将 瞬 态 RFI 识 别 为 子 
序列 ， 并 能 从 观测 数据 中 提取 瞬时 RFI 序 列 。 
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本 文 拟 提出 一 种 针对 瞬时 RFI 的 特征 提取 及 识别 方法 ， 期 望 针 对 特定 瞬时 RFI 实 现 特征 
提取 及 有 效 识别 。 


1 RF1 信 号 特征 提取 算法 


1.1 信号 周期 计算 

从 时 间 尺 度 上 观察 大 部 分 RFI 人 信号， 呈现 明显 的 周期 波动 。 为 了 分 析 RFI 信 号 的 特征 信 
Eh, 需要 对 信号 分 窗 处 理 ， 窗 口 大 小 与 信号 周期 信息 密切 相关 ,计算 信号 周期 是 分 析 其 特征 
信息 的 重要 前 提 。 

信号 周期 计算 算法 有 传统 的 倒 谱 法 、 自 相关 函数 法 、 平 均 幅 度 差 函数 法 和 小 波 变换 法 等 
PI。 自 相关 函数 法 是 一 种 非常 有 效 的 低 信 品 比 信号 周期 估计 算法 中 。 基 于 自 相 关 函 数 法 ， 
本 文 设计 了 一 种 自 相关 卷 积 滑 动 算法 计算 RFI 信 号 周期 。 算 法 过 程 如 下 : 

(1) ”对 一 段 长 度 为 N 的 RFI 信 号 振幅 值 取 绝对 值 。 

(2) ”从 RFI 信 号 上 截取 一 段 长 度 为 M 的 信号 作为 滑动 窗口 信号 。 

G) ”滑动 窗口 与 RFI 信 号 从 起 始 位 置 进行 卷 积 操作 ， 每 次 向 后 滑动 d 个 采样 点 。 

(4) ERP G) 直至 滑动 窗口 移动 到 RFI 信 号 尾部 。 

(5) ”利用 卷 积 后 结果 的 峰值 信息 ， 计 算 信 号 的 周期 大 小 。 


m 204 
ù | | 
a 
x 0 川 | 
Qa 
E | 
| | 
0.65 0.70 0.75 0.80 0.85 0.90 0.95 
Time(s) 
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Fig.1 (a) Time domain diagram of k7ing-3545khz noise signal 
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图 1 (b) k7ing-3545khz 噪 声 信号 周期 图 


Fig.1(b) Periodogram of k7ing-3545khz noise signal 


滑动 窗口 大 小 M 的 取 值 需要 小 于 RFI 信 号 长 度 和 NN ， 图 1 中 滑动 窗口 大 小 M -EN 为 
了 精确 地 计算 周期 ， 每 次 滑动 步 长 4d 取 值 为 1， 当 滑动 窗口 滑动 到 RFI 信 和 号 尾部 时 ， 剩 余 信 
号 长 度 如 果 小 于 滑动 窗口 信号 长 度 ， 需 要 进行 填 零 处 理 。 图 1 是 对 K7ing-3545khz' (ARRL 官 
网 mp3 格 式 RFI 信 号 ) 周期 分 机， 选取 了 如 图 1 (a) 0.65s~0.94s 约 0.29s 的 信号 。 周 期 计算 的 
结果 如 图 1 (b) 所 示 ， 可 以 明显 地 看 到 周期 性 尖峰 ， 通 过 求 平均 ， 计 算出 该 RFI 信 号 的 周期 。 
1. 2 RFI 信号 特征 提取 


1 http://www.arrl.org/files/file/RFI%20Sounds/k7ing-3545khz.mp3 
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是 RFI 信 号 重 


体 步 骤 如 下 : 
选取 RFI 
对 该 窗口 


要 的 特征 信息 ， 能 够 刻画 出 信号 波动 的 变化 细节 。RFI 信 号 特征 提取 


asain. 


a A) Se eA a EL 


获取 该 信号 窗口 的 峰值 信息 。 


对 


离散 峰值 信息 进行 预 处 理 ， 


使 其 平滑 。 


处 理 后 的 峰值 信息 


\ 作 为 RFI 信 号 的 特征 信息 。 


算法 


时 间 跨 度 


第 一 步 选取 RFI 信 号 
尽 可 能 短 且 至 少 包含 一 个 完整 RFI 周 期 信息 的 信 


周期 计生 
算法 将 信和 号 
周期 相关 ， 


法 确定 信号 周期 。RFI 信 号 
振幅 值 突 增 的 位 置 作为 窗口 的 起 始点 ，RFI 信 
计算 机 显示 器 -产生 的 干扰 信号 窗口 选择 如 图 2 所 示 ， 窗 
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窗口 ， 对 一 段 包含 RFI 的 信号 提取 特征 
号 段 ， 根 据 1.1 节 中 提出 的 RFI 信 和 号 
相对 于 其 他 信和 号， 振幅 值 较 大 ， 
FRIETEN 


信息 时 ， 需 要 选取 一 段 


起 始 和 终止 位 置 明确 ， 
取 窗 口 大 小 与 计算 得 到 的 
口 大 小 为 刀 。 
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图 2 RFI 信 号 窗口 选取 
Fig. 2 Selection of RFI signal window 
第 二 步 对 信号 的 振幅 值 求 绝对 值 ， 以 获取 更 多 的 峰值 信息 。 
第 三 步 利 用 振幅 值 特性 ， 获 取 窗 口 的 峰值 信息 。 
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Fig.3 (a) Peak curve of RFI signal 
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图 3 (b) 平滑 后 RFI 信 号 峰值 曲线 


Fig.3 (b) Smoothed RFI signal peak curve 


> http://www.arrl.org/files/file/RFI Sounds/monitor.mp3 
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第 四 步 对 获取 的 峰值 信号 进行 预 处 理 , A IREK Ae Fal BA EY EL ES F 


滑 峰值 曲线 ， 减 少 局 部 波动 ， 纵 小 同类 别 信号 差异 ， 提 高 识别 效率 。 图 3 a 为 RFI 信 号 峰 


值 信息 折线 图 , 波动 趋势 不 平稳 ， 不 利于 后 续 特 征 识别 。 进 行 平滑 处 理 后 ， 如 图 3(b) 所 示 ， 
曲线 相对 平稳 ， 能 够 缩小 同类 信号 的 差异 。 将 平滑 处 理 后 的 峰值 信息 作为 识别 RFI 信 号 的 特 


征 模 板 ， 可 进行 相似 度 计算 。 


2 RF1 信 号 特征 相似 度 算法 


对 未 知 信号 进行 特征 提取 后 ， 计 算 与 已 知 特征 模板 的 相似 程度 ， 实 现 对 未 知 RFI 信 号 的 


识别 和 分 类 。 算 法 采用 分 段 打 分 策略 


可 以 分 为 多 个 区 段 ， 比 较 两 个 序列 对 应 区 段 的 趋势 是 否 相 同 ， 进 行 记 分 ,趋势 相 同 加 分 ， 否 


则 减 分 。 

在 进行 相似 度 计 算 前 ， 需 要 确保 
否 一 一 对 应 。 竺 识别 信号 具体 细节 未 
征 模板 对 齐 ， 即 相位 信息 没有 对 章 ， 
做 准备 。 


计算 两 个 序列 的 相似 程度 , 每 一 个 序列 根据 相 邻 离散 点 


这 两 个 序列 处 于 “对 齐 ” 状 态 ， 即 确定 两 个 序列 的 区 段 是 
Al, 提取 峰 值 信息 后 ， 离 散 的 峰值 信息 几乎 不 可 能 与 特 
因此 其 峰值 位 置信 息 等 都 需要 重新 调整 ， 为 计算 相似 度 


文献 [13] 提 出 了 一 种 测量 时 间 序 列 相似 性 的 方法 ， 即 动态 时 间 规 整 《Dynamic Time 


Warping, DTW) ， 该 方法 能 够 比较 


不 同 长 度 的 时 间 序列 0， 已 广泛 应 用 于 语音 信号 处 理 ， 


在 该 算法 中 计算 得 到 的 欧 几 里 得 距离 越 小 , 表明 两 种 声音 模式 的 相似 度 越 高 " 。 利用 DTW 


算法 能 够 计算 出 两 个 长 度 不 同 序列 的 欧 几 里 得 距离 最 小 时 , 这 两 个 序列 离散 点 对 应 情况 ,如 
图 4 所 示 ， 两 个 序列 的 相位 点 对 应 情况 。 


Fig. 4 Correspondence between the phas 


图 4 使 用 DTW 算 法 后 两 个 序列 的 相位 点 对 应 情况 


e points of the two sequences after using the DTW algorithm 


基于 DTW 算法 ， 本 文 设计 了 RFI 信号 特征 相似 度 计 算 算法 ， 其 过 程 如 下 : 


(1) ”通过 DTW 算法 计算 两 组 序列 在 最 短 欧 几 里 得 距离 条 件 下 ， 离 散 点 对 应 情况 。 
(2) ”计算 未 知 信 号 离散 点 的 权重 值 。 
G) ”分 段 比 较 ， 如 果 对 应 段 的 变化 趋势 相同 ， 相 似 度 分 数 增加 ， 如 果 变 化 趋势 不 同 ， 


则 相似 度 分 数 减少 。 


(4) ”累计 各 段 的 分 数 ， 总 分 数 代表 两 个 序列 的 相似 程度 。 


RFI 信 号 特征 相似 度 计算 算法 第 


离 下 离散 点 对 应 情况 ， 为 后 续 计 算 做 
如 公式 (1) 所 示 。 


步 利 用 DTW 算 法 求 得 待 比 较 序列 在 最 小 欧 几 里 得 距 
准备 。 第 二 步 计算 被 比较 序列 的 离散 点 权 值 ， 计 算 方 法 


3 https://en.wikipedia.org/wiki/Euclidean_distance 


SO ogg 
WOOD=1> SH (GD 
0 i=l or i=D 
HE DAR SIMA BUR REL S(t) 代表 该 序列 第 t 个 离散 点 的 幅 值 。W(7) 是 序列 第 t 个 


离散 点 的 权重 值 ， 离 散 值 越 大 ， 权 重 越 大 ,在 信号 特征 中 的 代表 性 越 强 。 为 了 使 相似 度 的 大 
小 介 于 -1 到 1 之 间 ， 第 一 个 和 最 后 一 个 的 权重 设置 为 0。 

第 三 步 进行 序列 相似 度 计算 ， 采 用 分 段 打分 累计 方式 。 通 过 DTW 算 法 计算 每 一 区 段 的 
对 应 情况 ， 采 用 分 段 比 较 方 式 ， 如 图 4 所 示 两 个 序列 的 区 段 有 趋势 相同 、 趋 势 不 同 和 一 对 多 
三 种 情况 ， 两 个 区 段 的 趋势 相同 时 加 分 ， 如 公式 〈2) Aras. 


score = score + M [i] (2) 


其 中 score 是 两 个 序列 相似 度 的 累积 分 数 。AMX 四 是 第 ;个 区 段 的 权重 ， 值 为 区 段 的 两 个 离散 


点 权重 的 平均 值 ， 即 0.5x《W[]+W[i+1]) 。 两 个 区 段 的 趋势 不 同时 减 分 ， 如 公式 〈3) 所 


ZN o 


score = score — M [i] (3) 


对 于 一 对 多 情况 不 进行 加 减 分。 计算 并 累计 各 段 的 分 数 后 得 到 两 个 序列 相似 度 值 。 


3 实验 结果 与 分 析 


实验 选择 ARRL 官 网 "中 提供 的 干扰 源 作为 实验 数据 ， 选 择 七 种 REFI 源 进行 交叉 测试 ， 基 
本 信息 如 表 1 所 示 。 为 方便 对 RFI 信 号 进行 特征 提取 及 后 续 识 别 测试 ， 利 用 1.1 节 提出 的 RFI 
童 号 周期 计算 法 得 到 各 实验 数据 周期 采样 点 数 ， 以 便 下 一 步 特征 提取 及 识别 测试 。 


表 1 RFI 源 的 基本 信息 


Tab. 1 Basic information on RFI sources 


REFI 来 源 采样 率 (Hz) ”每 周期 采样 点 数 ”周期 (s) 
k7ing-3545khz 32000 529 0.01653 
ks2am-streetlight 32000 547 0.01709 
monitor 22050 372 0.01687 
n6rce-sps-carrier 44100 728 0.01650 
ple-4 32000 1321 0.04128 
ausoth 44100 797 0.01807 
181200th25 11025 443 0.04018 


4 http://www.arrl.org/sounds-of-rfi 
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随机 选取 每 个 信号 的 10 个 周期 作为 测试 数据 ,提取 出 相应 候选 特征 模板 与 其 他 信号 进行 
交叉 测试 , 选择 与 相同 信号 相似 度 大 且 与 其 他 信号 相似 度 较 低 的 候选 模板 为 最 优 模板 。 如 图 
5 (a) 选取 K7_2 作 为 K7ing-3545khz 的 特征 模板 ， 图 3 (b) 选取 ks_9 作 为 ks2am-streetlight 的 特 
征 模 板 ， 图 $〈c) 选取 mo_4 作 为 monitor 的 特征 模板 ， 图 5 Cd) 选取 n6_10 作 为 n6rce-sps-carrier 
的 特征 模板 ， 图 5(e) 选取 pl_8 作 为 ple-4 的 特征 模板 ， 图 5 Cf) 选取 au_7 作 为 ausoth 的 特征 
Bi, PAS Ce) 选取 ot 4 作为 18120oth25 的 特征 模板 。 
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图 5 Ca) 选择 k7ing-3545khz 不 同 模板 的 相似 度 结果 
Fig.5 (a) Similarity results for different templates selected for k7ing-3545khz 
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图 5 (b) 选择 ks2am-streetlight 不 同 模板 的 相似 度 结果 


Fig.5 (b) Similarity results for different templates selected for ks2am-streetlight 
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Scores by group 
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Fig.5 (c) Similarity results for different templates selected for monitor 
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图 5 Cd) 选择 n6rce-sps-carrier 不 同 模板 的 相似 度 结果 


Fig.5 (d) Similarity results for different templates selected for n6rce-sps-carrier 
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图 5〈e) 选择 plc-4 不 同 模板 的 相似 度 结果 


Fig.5 Ce) Similarity results for different templates selected for plc-4 
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Scores by group 
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图 5 (f) 选择 ausoth 不 同 模板 的 相似 度 结果 


Fig.5 (f) Similarity results for different templates selected for ausoth 
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KIS Cg) 选择 18120oth25 不 同 模板 的 相似 度 结果 


Fig.5 (g) Similarity results for different templates selected for18120oth25 


表 2 RFI 信 号 相似 度 验 证 平均 值 结 果 


Tab.2 RFI signal similarity verification average results 


特征 模板 k7ing-3545khz ks2am-streetlight monitor n6rce-sps-carrier ple-4 ausoth 181200th25 
k7ing-3545khz 0.875 0.107 0.150 0.133 0.335 0.073 0.142 
ks2am-streetlight 0.263 0.723 0.549 0.277 0.446 0.582 0.403 
monitor 0.468 0.583 0.792 0.478 0.566 0.520 0.637 
n6rce-sps-catrier 0.445 0.479 0.571 0.893 0.726 0.400 0.449 
plc-4 0.412 0.172 0.205 0.169 0.762 0.057 0.338 
ausoth 0.187 0.479 0.395 0.180 0.211 0.870 0.488 
181200th25 0.262 0.299 0.344 0.178 0.522 0.270 0.881 


提取 每 个 信号 最 优 特征 模板 所 在 组 数据 ， 结 果 如 表 2 所 示 ， 相 同 信号 的 相似 度 在 所 在 行 
最 大 ， 表 明 提取 出 的 RFI 特 征 模板 能 够 较 正 确 地 识别 RFI 来 源 。 算 法 的 相似 度 计 算 基 于 特征 
模板 得 出 ， 权 值 计 算 依赖 于 特征 模板 ， 即 计算 两 个 序列 的 相似 度 时 ， 权 值 的 大 小 只 与 特征 模 
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板 有 关 ， 因 此 以 对 角 线 互 为 对 称 的 数值 不 同 。 同 一 信号 在 不 同 周期 内 存在 差别 ， 造 成 同一 信 
号 的 相似 度 小 于 1。 

表 2 中 个 别 数 据 偏 高 ， 如 n6rce-sps-carrier 行 plc-4 列 的 数值 为 0.726， 通 过 分 析 得 知 待 比较 
的 一 组 序列 振幅 波动 趋势 比较 一 致 ， 导 致 相似 度 偏 高 ， 后 续 的 工作 将 继续 优化 算法 ， 以 提高 
其 精度 。 


4 结语 


本 文 设计 了 滑动 卷 积 周期 计算 算法 ， 利 用 卷 积 后 的 峰值 区 间 完 成 了 信和 号 平均 周期 计算 ; 
基于 设计 的 峰值 提取 算法 实现 了 信和 号 特征 模板 提取 ;， 基于 DTW 算 法 和 打分 策略 设计 了 信和 号 
特征 识别 算法 ， 实 现 了 未 知 信号 的 识别 和 分 类 。 分 别 对 ARRL 官 网 下 载 的 REFI 数 据 进行 了 互 
相关 与 自 相关 计算 , 实验 结果 表明 , 信号 生成 的 模板 与 原始 信号 的 相似 度 明 显 高 于 其 他 信号， 
说 明 本 文 提出 的 算法 可 以 有 效 地 提取 特征 信息 并 生成 信号 特征 模板 ， 且 利用 RFI 特 征 识别 算 
法 能 对 RFI 进 行 正确 分 类 。 本 文 提 出 的 方法 可 对 RFI 进 行 细 粒度 识别 和 标记 ， 有 望 为 射电 天 
文 RFI 特 征 识别 和 标记 提供 一 种 新 的 解决 方案 。 

本 论文 中 算法 实现 代码 以 及 实验 数据 已 在 码 云 仓库 开源 。 
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Abstract: Transient radio frequency interference (RFI) signals are one of the most difficult 
types of interference to identify and eliminate. This paper proposes a feature extraction method 
based on RFI signal peaks for the problem of transient RFI encountered during radio astronomical 
observations. The method is used to generate a feature template for RFI signals, and a feature 
recognition algorithm is designed based on the Dynamic Time Warping (DTW) algorithm and 
scoring strategy to achieve RFI signal recognition and classification. The cross-validation 
comparison of the feature templates was completed using the RFI data provided by the ARRL 
website. The test results show that the RFI feature templates are highly similar to the same 
category of RFI, indicating that the algorithm is capable of identifying and classifying RFI, which 
provides a new idea for the feature extraction and identification of transient RFI. 
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